Nematic topological defects in the presence of axisymmetric fluid flow 
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Recent numerical simulations of lid-driven cavity flow of a nematic liquid crystal have found 
dynamical behavior where topological defects rotate about the center of the fluid vortex induced by 
the lid motion. By considering a simpler geometry of an infinite system with axisymmetric fluid flow 
we show that the Ericksen-Leslie nematodynamic equation for the director can be solved exactly. 
The solution demonstrates that any configuration of defects will be advected by the fluid flow, with 
the defects rotating about the center of the fluid vortex with the angular velocity of the fluid. 

PACS numbers: 61.30Dk,61.30Jf,61.30Gd 



Lid-driven cavity flow is a standard geometry used 
in both experimental and computational studies of flu- 
ids. In this geometry a fluid is contained within a box 
whose sides are stationary except for the the top wall 
(the "lid") which is driven at some nonzero velocity. In 
numerical studies the problem is often simplified to a 
two-dimensional square domain with one moving side. 
The driven lid creates a vortex fluid flow near the center 
of the cavity as well as counterrotating vortices near the 
corners. If the fluid is a complex one, such as a nematic 
liquid crystal, the driven lid can lead to very interest- 
ing and rich topological defect phenomena as has been 
observed recently in a number of numerical simulation 
studies [l]-[H|- Among the phenomena observed are ro- 
tating half-integer charged topological defects, either a 
single defect or a pair of oppositely charged defects 
[l[ , depending on the director boundary conditions em- 
ployed. The defects are advected by the fluid vortex near 
the center of the cavity and undergo a periodic rotation 
with the fluid flow about the center of the vortex. In the 
case of a single defect Q the defect core eventually settles 
at the center of the vortex and remains stationary, while 
in the case of a defect pair [l[ the defects continue to ro- 
tate indefinitely with other transient defect phenomena 
occuring elsewhere in the cavity. 

In this paper we show that by considering a sim- 
pler but related geometry the Ericksen-Leslie nemato- 
dynamic equation for the director can be solved exactly 
for any texture composed of an array of defects. The 
solution is consistent with the behavior observed in the 
simulations, namely, advection of the defects about the 
center of the fluid vortex with an angular speed equal to 
that of the fluid. 

Specifically, we consider a nematic liquid crystal with 
a simple vortex fluid flow (a "stationary axisymmetric 
flow" ) given in plane polar coordinates by: 



(1) 



where r measures the radial distance from the center of 
the vortex and is the tangential polar unit vector. 
The angular speed of the fluid ui is the magnitude of the 
angular velocity: 

oj = ly x v. (2) 



We assume that the system is infinite in size, thus sim- 
plifying the cavity flow geometry while retaining the 
essence of the lid-driven flow, namely, the central vor- 
tex. For this vortex flow the symmetric velocity tensor 

is identically zero and the Ericksen- 
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Leslie hydrodynamic equation for the nematic director 
simplifies to: 

7n = K\7 2 n + (a 3 - a 2 )N 



(3) 



assuming a single Frank elastic constant K. Here ot2,U3 
are Leslie viscosity coefficients and N is the rotation rate 
of the director relative to the moving fluid: 
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v • Vn — lj x n 



(4) 



The Lagrange multiplier 7 maintains the constraint: 
|n| = l. 

We assume that the director field lies in the plane of the 
fluid flow (the x — y plane) and express it in rectangular 
coordinates by: 



(cos9(x, y), sin9(x, y)) 



(5) 



where 9 is the angle n makes with the x axis. Substitut- 
ing Eqs. (JTJ) and ([5]) into Eq. ([3]) the two nonzero compo- 
nents of the latter equation are found to be: 
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where A = «3 — oi2- After eliminating the unknown La- 
grange multiplier 7 from the above two equations we ob- 
tain: 

a* + «(*£-'£)'-"' -3™ (9) 

The left hand side of Eq. © is the time rate of change 
of the director angle 9 relative to the rotating fluid [j| . 

If initially V 2 6* = 0, a condition satisfied at all points in 
space by any configuration of topological defects except 
at the isolated defect cores 0, then Eq. © requires that 
the time rate of change of 8 relative to the rotating fluid 
is zero, i.e., 9 as a function of the polar coordinates in 
the nonrotating frame must satisfy: 

9(r, cj>, t) = 9(r, (j> - uAt, t - At) + uAt (10) 

for all t and At Thus, if there are topological defects in 
the system initially the director pattern will rotate with 
the fluid, i.e., the defects will rotate with same angular 
velocity uj as the fluid. The above equation is the exact 
solution to the Ericksen-Leslie director equation for an 
axisymmetric flow and an initial director configuration 
satisfying V 2 9 — 0. While the latter equation is not 
satisfied at the defect core itself, the core must rotate at 
the angular velocity as the director field surrounding it 
to preserve the continuity of the structure. Also, note 
that this solution is stable against director fluctuations 
out of the plane as can be seen from looking at the z 
component of Eqn. ([3]). 



In Figs. Q] and [2] we illustrate Eq. (IT0|) for two cases: a 
single half-integer defect and a pair of oppositely charged 
half-integer defects, configurations that were observed 
in the numerical simulations of lid-driven cavity flow in 
Refs. @ and [l[, respectively. In both figures we have 
chosen uj = — 7r, where the negative sign corresponds to 
the clockwise fluid flow imposed in the simulations by 
a lid driven in the positive x direction. The figures il- 
lustrate the periodic rotation of the defects with period 
2n /lj = 2 as observed in the numerical simulations. This 
is the steady-state motion for the axisymmetric infinite 
system considered in this paper. In the simulation of 
Ref. [|[ the single defect eventually settled to the center 
of the vortex; this behavior arises presumably from the 
presence of the defects at the corners of the cavity leading 
to an equilibrium position of the rotating defect at the 
vortex center. In the case of the rotating pair of defects 
1] the rotation was observed to occur indefinitely as is 
the case for our simple geometry. 
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FIG. 1: Illustration of the motion of a single half-integer topological defect with the director angle 6 given by Eq. (|10[) . We 
have chosen oj = — it, i.e, a clockwise fluid flow and placed the defect initially at x = —4.3. Compare with Fig. 5 of [J 
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(e) t = 1.6 (f) t = 2.0 

FIG. 2: Illustration of the motion of a pair of oppositely charged half-integer topological defects with the director angle 9 given 
by Eq. (|10[) . We have chosen uj = —n, i.e, a clockwise fluid flow and placed the defects initially at a distance 13.5 on either side 
of the origin. Compare with Fig. 8 of [l|. 



